Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

The de novo genome of the Black-necked Snakefly (Venustoraphidia nigricollis Albarda, 1891): A resource to study the evolution of living fossils.

Snakeflies (Raphidioptera) are the smallest order of holometabolous insects that have kept their distinct and name-giving appearance since the Mesozoic, probably since the Jurassic, and possibly even since their emergence in the Carboniferous, more than 300 million years ago. Despite their interesting nature and numerous publications on their morphology, taxonomy, systematics, and biogeography, snakeflies have never received much attention from the general public, and only a few studies were devoted to their molecular biology. Due to this lack of molecular data, it is therefore unknown, if the conserved morphological nature of these living fossils translates to conserved genomic structures. Here, we present the first genome of the species and of the entire order of Raphidioptera. The final genome assembly has a total length of 669 Mbp and reached a high continuity with an N50 of 5.07 Mbp. Further quality controls also indicate a high completeness and no meaningful contamination. The newly generated data was used in a large-scaled phylogenetic analysis of snakeflies using shared orthologous sequences. Quartet score and gene concordance analyses revealed high amounts of conflicting signals within this group that might speak for substantial incomplete lineage sorting and introgression after their presumed re-radiation after the asteroid impact 66 million years ago. Overall, this reference genome will be a door-opening dataset for many future research applications, and we demonstrated its utility in a phylogenetic analysis that provides new insights into the evolution of this group of living fossils.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality.

BACKGROUND: Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS: Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS: Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae).

BACKGROUND: Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics. RESULTS: The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function. CONCLUSIONS: Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

A chromosome-scale high-contiguity genome assembly of the cheetah (Acinonyx jubatus).

The cheetah (Acinonyx jubatus, SCHREBER 1775) is a large felid and is considered the fastest land animal. Historically, it inhabited open grassland across Africa, the Arabian Peninsula, and southwestern Asia; however, only small and fragmented populations remain today. Here, we present a de novo genome assembly of the cheetah based on PacBio continuous long reads and Hi-C proximity ligation data. The final assembly (VMU_Ajub_asm_v1.0) has a total length of 2.38 Gb, of which 99.7% are anchored into the expected 19 chromosome-scale scaffolds. The contig and scaffold N50 values of 96.8 Mb and 144.4 Mb, respectively, a BUSCO completeness of 95.4% and a k-mer completeness of 98.4%, emphasize the high quality of the assembly. Furthermore, annotation of the assembly identified 23,622 genes and a repeat content of 40.4%. This new highly contiguous and chromosome-scale assembly will greatly benefit conservation and evolutionary genomic analyses and will be a valuable resource, e.g., to gain a detailed understanding of the function and diversity of immune response genes in felids.

The potential and shortcomings of mitochondrial DNA analysis for cheetah conservation management.

There are only about 7,100 adolescent and adult cheetahs (Acinonyx jubatus) remaining in the wild. With the majority occurring outside protected areas, their numbers are rapidly declining. Evidence-based conservation measures are essential for the survival of this species. Genetic data is routinely used to inform conservation strategies, e.g., by establishing conservation units (CU). A commonly used marker in conservation genetics is mitochondrial DNA (mtDNA). Here, we investigated the cheetah's phylogeography using a large-scale mtDNA data set to refine subspecies distributions and better assign individuals to CUs. Our dataset mostly consisted of historic samples to cover the cheetah's whole range as the species has been extinct in most of its former distribution. While our genetic data largely agree with geography-based subspecies assignments, several geographic regions show conflicting mtDNA signals. Our analyses support previous findings that evolutionary forces such as incomplete lineage sorting or mitochondrial capture likely confound the mitochondrial phylogeography of this species, especially in East and, to some extent, in Northeast Africa. We caution that subspecies assignments solely based on mtDNA should be treated carefully and argue for an additional standardized nuclear single nucleotide polymorphism (SNP) marker set for subspecies identification and monitoring. However, the detection of the A. j. soemmeringii specific haplogroup by a newly designed Amplification-Refractory Mutation System (ARMS) can already provide support for conservation measures. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01483-1.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding.

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics.

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

Virus diversity in metagenomes of a lichen symbiosis (Umbilicaria phaea): complete viral genomes, putative hosts and elevational distributions.

Viruses can play critical roles in symbioses by initiating horizontal gene transfer, affecting host phenotypes, or expanding their host's ecological niche. However, knowledge of viral diversity and distribution in symbiotic organisms remains elusive. Here we use deep-sequenced metagenomic DNA (PacBio Sequel II; two individuals), paired with a population genomics approach (Pool-seq; 11 populations, 550 individuals) to understand viral distributions in the lichen Umbilicaria phaea. We assess (i) viral diversity in lichen thalli, (ii) putative viral hosts (fungi, algae, bacteria) and (iii) viral distributions along two replicated elevation gradients. We identified five novel viruses, showing 28%-40% amino acid identity to known viruses. They tentatively belong to the families Caulimoviridae, Myoviridae, Podoviridae and Siphoviridae. Our analysis suggests that the Caulimovirus is associated with green algal photobionts (Trebouxia) of the lichen, and the remaining viruses with bacterial hosts. We did not detect viral sequences in the mycobiont. Caulimovirus abundance decreased with increasing elevation, a pattern reflected by a specific algal lineage hosting this virus. Bacteriophages showed population-specific patterns. Our work provides the first comprehensive insights into viruses associated with a lichen holobiont and suggests an interplay of viral hosts and environment in structuring viral distributions.

Two high-quality de novo genomes from single ethanol-preserved specimens of tiny metazoans (Collembola).

BACKGROUND: Genome sequencing of all known eukaryotes on Earth promises unprecedented advances in biological sciences and in biodiversity-related applied fields such as environmental management and natural product research. Advances in long-read DNA sequencing make it feasible to generate high-quality genomes for many non-genetic model species. However, long-read sequencing today relies on sizable quantities of high-quality, high molecular weight DNA, which is mostly obtained from fresh tissues. This is a challenge for biodiversity genomics of most metazoan species, which are tiny and need to be preserved immediately after collection. Here we present de novo genomes of 2 species of submillimeter Collembola. For each, we prepared the sequencing library from high molecular weight DNA extracted from a single specimen and using a novel ultra-low input protocol from Pacific Biosciences. This protocol requires a DNA input of only 5 ng, permitted by a whole-genome amplification step. RESULTS: The 2 assembled genomes have N50 values >5.5 and 8.5 Mb, respectively, and both contain ∼96% of BUSCO genes. Thus, they are highly contiguous and complete. The genomes are supported by an integrative taxonomy approach including placement in a genome-based phylogeny of Collembola and designation of a neotype for 1 of the species. Higher heterozygosity values are recorded in the more mobile species. Both species are devoid of the biosynthetic pathway for ß-lactam antibiotics known in several Collembola, confirming the tight correlation of antibiotic synthesis with the species way of life. CONCLUSIONS: It is now possible to generate high-quality genomes from single specimens of minute, field-preserved metazoans, exceeding the minimum contig N50 (1 Mb) required by the Earth BioGenome Project.

Inference of DNA methylation patterns in molluscs.

Mollusca are the second largest and arguably most diverse phylum of the animal kingdom. This is in sharp contrast to our very limited knowledge concerning epigenetic mechanisms including DNA methylation in this invertebrate group. Here, we inferred DNA methylation patterns by analysing the normalized dinucleotide CG content in protein-coding sequences and identified DNA methyltransferases (DNMT1 and 3) in published transcriptomes and genomes of 140 species across all eight classes of molluscs. Given the evolutionary age and morphological diversity of molluscs, we expected to find evidence for diverse methylation patterns. Our inferences suggest that molluscs possess substantial levels of DNA methylation in gene bodies as a rule. Yet, we found deviations from this general picture with regard to (i) the CpG observed/expected distributions indicating a reduction in DNA methylation in certain groups and (ii) the completeness of the DNMT toolkit. Reductions were evident in Caudofoveata, Solenogastres, Polyplacophora, Monoplacophora, as well as Scaphopoda. Heterobranchia and Oegopsida were remarkable as they lacked DNMT3, usually responsible for de novo methylation, yet showed signs of DNA methylation. Our survey may serve as guidance for direct empirical analyses of DNA methylation in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

A new genus and species of mud snake from Myanmar (Reptilia, Squamata, Homalopsidae).

Based on two male and two female individuals, we describe a new genus and species of mud snake, Myanophis thanlyinensis gen. nov., sp. nov., from the vicinity of the campus of East Yangon University, Yangon, Thanlyin, Myanmar. This species differs from every other homalopsid species by the following combination of characters: (1) dorsal scales smooth, row formula 21-21-19 or 21-21-17; (2) tail short, ratio tail length/SVL 0.185-0.204 in males, 0.160-0.167 in females; (3) nasal scales separated; (4) 125-126 ventral scales in males, 120-122 in females; (5) 38-39 subcaudal scales in males, 32-34 in females; and (6) hemipenis bilobed. Its matrilineal genealogy (based on analyses of 16S and cytochrome b sequences), associates Myanophis thanlyinensis gen. nov., sp. nov. most closely with species of the genera Myrrophis and Gyiophis. The new taxon differs from the species of Myrrophis and Gyiophis by having a bilobed hemipenis (vs. unilobed). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Myrrophis by having 125-126 ventral scales in males and 120-122 in females (vs. 137-162 and 137-164, respectively), and 38-39 subcaudal scales in males and 32-34 in females (vs. 39-55 and 37-52, respectively). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Gyiophis by lacking dark blotches along flank (vs. present), and by having 21 dorsal scales rows at midbody (vs. 25). We provide an identification key to the homalopsid species known to occur in Myanmar. As a novelty to the classic holotype description and characterization, the individual has been genome sequenced by Illumina short-read technology and its genome has been assembled into a draft nuclear genome and a complete, annotated mitochondrial genome. This innovative approach comprehensively and permanently characterizes the genomic variation of the holotype.

Detection of SARS-CoV-2 in raw and treated wastewater in Germany - Suitability for COVID-19 surveillance and potential transmission risks.

Wastewater-based monitoring of the spread of the new SARS-CoV-2 virus, also referred to as wastewater-based epidemiology (WBE), has been suggested as a tool to support epidemiology. An extensive sampling campaign, including nine municipal wastewater treatment plants, has been conducted in different cities of the Federal State of North Rhine-Westphalia (Germany) on the same day in April 2020, close to the first peak of the corona crisis. Samples were processed and analysed for a set of SARS-CoV-2-specific genes, as well as pan-genotypic gene sequences also covering other coronavirus types, using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Additionally, a comprehensive set of chemical reference parameters and bioindicators was analysed to characterize the wastewater quality and composition. Results of the RT-qPCR based gene analysis indicate the presence of SARS-CoV-2 genetic traces in different raw wastewaters. Furthermore, selected samples have been sequenced using Sanger technology to confirm the specificity of the RT-qPCR and the origin of the coronavirus. A comparison of the particle-bound and the dissolved portion of SARS-CoV-2 virus genes shows that quantifications must not neglect the solid-phase reservoir. The infectivity of the raw wastewater has also been assessed by viral outgrowth assay with a potential SARS-CoV-2 host cell line in vitro, which were not infected when exposed to the samples. This first evidence suggests that wastewater might be no major route for transmission to humans. Our findings draw attention to the need for further methodological and molecular assay validation for enveloped viruses in wastewater.

Correction to: An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola).

An amendment to this paper has been published and can be accessed via the original article.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola).

BACKGROUND: The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations. RESULTS: Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera. CONCLUSION: Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA.

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

Snails in the desert: Species diversification of Theba (Gastropoda: Helicidae) along the Atlantic coast of NW Africa.

The spatial subdivision of species often plays a pivotal role in speciation. Across their entire range, species are rarely panmictic and crucial consequences of spatial subdivision are (1) random genetic drift including historical factors, (2) uniform selection, and (3) divergent selection. Each of these consequences may result in geographic variation and eventually reproductive isolation, but their relative importance in speciation is still unclear. In this study, we used a combination of genetic, morphological, and climatic data to obtain a comprehensive picture of differentiation among three closely related, parapatrically distributed taxa of the land snail genus Theba occurring along the Atlantic coasts of South Morocco and Western Sahara. We conducted Mantel and partial Mantel tests to relate phenotypic and genotypic variation of these species to geography and/or climate. As null hypothesis for an evolutionary scenario, we assumed nonadaptive speciation and expected a pattern of isolation by distance among taxa. Rejection of the null hypothesis would indicate isolation by environment due to adaptation. Generally, genetic drift plays an important role but is rarely considered as sole driver of speciation. It is the combination of drift and selection that predominantly drives speciation. This study, however, provides a potential example, in which nonadaptive speciation, that is, genetic drift, is apparently the main driver of shaping the diversity of Theba in NW Africa. Restriction of gene flow between populations caused by geographic isolation probably has played an important role. Climate oscillations during the Plio- and Pleistocene may have led to repeated ecological changes in NW Africa and disruptions of habitats promoting differentiation by geographic isolation. The inferred evolutionary scenario, however, did not fully explain the incongruence between the AFLP- and mtDNA-tree topologies. This incongruence might indicate past hybridization among the studied Theba forms.

The complete mitochondrial genome of the 'solar-powered' sea slug Plakobranchus cf. ocellatus (Heterobranchia: Panpulmonata: Sacoglossa).

We present the complete mitochondrial genome sequence of Plakobranchus cf. ocellatus (Heterobranchia: Sacoglossa), a so-called 'solar-powered' sea slug with long-term retention of chloroplasts. The mitochondrial genome was 14,177 bp in length containing the standard set of 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The base composition of 27.3% A, 15.6% C, 18.6% G, and 38.5% T showed a strong A + T bias. The genome organization of P. cf. ocellatus is identical to the other sacoglossan mitogenomes sequenced so far, except for Ascobulla fragilis.

Phylogenomics resolves the timing and pattern of insect evolution.

Insects are the most speciose group of animals, but the phylogenetic relationships of many major lineages remain unresolved. We inferred the phylogeny of insects from 1478 protein-coding genes. Phylogenomic analyses of nucleotide and amino acid sequences, with site-specific nucleotide or domain-specific amino acid substitution models, produced statistically robust and congruent results resolving previously controversial phylogenetic relations hips. We dated the origin of insects to the Early Ordovician [~479 million years ago (Ma)], of insect flight to the Early Devonian (~406 Ma), of major extant lineages to the Mississippian (~345 Ma), and the major diversification of holometabolous insects to the Early Cretaceous. Our phylogenomic study provides a comprehensive reliable scaffold for future comparative analyses of evolutionary innovations among insects.